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Abstract 

We derive the phase space particle density operator in the 'droplet' 
picture of bosonization in terms of the boundary operator. We demon- 
strate that it satisfies the correct algebra and acts on the proper Hilbert 
space describing the underlying fermion system, and therefore it can be 
used to bosonize any hamiltonian or related operator. As a demonstration 
we show that it reproduces the correct excitation energies for a system of 
free fermions with arbitrary dispersion relations. 



1 Introduction 



The description of a one-dimensional fermion system in terms of bosonic degrees 
of freedom, called bosonization, is a well-established and useful technique ^ |2J 
E| • The advantage of such a transmutation lies in the fact that the dynamics 
of the equivalent bosonic theory encodes many-body excitations of the fermion 
system, with large coherent (classical) bosonic excitations representing large 
collective states of the fermions. It provides, therefore, a convenient setting to 
study strongly correlated fermions. 

It is intuitively obvious that such bosonic descriptions must trivially exist 
in any number of dimensions, at least within each fermion number sector: any 
excitation of the system that does not change the number of fermions (or changes 
them by an even number) is inherently bosonic and can be obtained with the 
action of a bosonic operator on the ground state. The full theory, therefore, 
can be described in terms of a collection of these bosonic operators. Examples 
of such operators are the "hydrodynamic" current and particle densities, or the 
more general phase space density operator that will be examined in this paper. 

Such descriptions, however, do not qualify as proper bosonizations. The set 
of bosonic operators describing the system would form, in general, an overcom- 
plete set and would make the identification of the proper features of the fermion 
system difficult. More formally, the algebra of the bosonic operators would, in 
general, admit several irreducible representations, of which only one would be 
the fermion Hilbert space. Consequently, the classical action for the bosonic 
fields could not be used to identify proper collective excitations of the fermions. 

The above-mentioned hydrodynamic or phase space density variables provide 
examples of such a redundancy: their algebra is the same for all systems, and 
admits representations that describe fermions, bosons, as well as a host of other 
systems (parastatistics, fractional filling states etc.) A successful bosonization 
of a system should realize the fermion states as its complete, irreducible Hilbert 
space. 

Bosonization works successfully for one spatial dimension, where special 
properties of line- like spaces come into play. In particular, statistics and spin are 
essentially irrelevant notions in one dimension and thus fermionic and bosonic 
theories could a priori have the same degrees of freedom. On a more intuitive 
level, fermion dynamics around the Fermi (or Dirac) sea become tantamount 
to one-dimensional wave propagation, with the corresponding phonon states 
encoding quantum excitations pQ. Further, exponentials in the bosonic field 
('vertex operators') can be properly defined to obey fermionic commutation re- 
lations. Overall, we have an exact mapping of states between the two systems 
both at the many-body and the field theory level. 

Bosonization can be motivated or formulated from various starting points 
(operator, field theory, many-body etc.). A particularly intuitive approach is 
to consider a semiclassical description in terms of a dense collection of fermions 
forming a 'droplet' in their phase space and study the dynamics of the Fermi 
surface of the droplet 0] . In one dimension, this leads to a chiral theory that fully 
reproduces the standard bosonization results. In particular, for nonrelativistic 
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fermions, it reproduces the results of the collective field method [7]. 

The basic advantage of the phase space approach is that, in principle, it 
works in any number of dimensions and it seems to capture the relevant degrees 
of freedom for a successful bosonization. In fact, an adaptation of this method 
properly taking into account the quantum nature of the phase space was pro- 
posed by one of the authors as the starting point for an exact (that is, complete 
and irreducible) bosonization in arbitrary dimensions [5]- This could become a 
useful tool in the study of correlated fermion systems in dimensions higher that 
one. (For other approaches to higher dimensional bosonization see 

In paving the way to such applications, it is important to establish the 
full description of the Fermi system in the phase space droplet picture. In 
particular, it is important to have an exact mapping of all fermion operators in 
terms of the (bosonic) field describing the droplet boundary. This would provide 
the bosonized form of any many-body hamiltonian, including interactions, and 
requires deriving the boundary field expression of a complete set of fermion 
operators. Such a set is encoded by the phase space density operator, in terms 
of which any other operator can be obtained. 

The purpose of this paper is to do this in the controllable case of one spatial 
dimension, in which standard bosonization works. Specifically, the expression 
for the quantum density operator will be given as a function of the droplet 
boundary field. It will be shown that this operator obeys the correct algebra 
and acts on the correct Hilbert space provided by the unique irreducible rep- 
resentation of the boundary field. It therefore realizes a bosonization of any 
fermion system in the droplet picture and can be used as the starting point for 
higher-dimensional generalizations . 

2 Review of the droplet method 

The starting point is a semiclassical description of the many-body fermion state 
in terms of its phase space density. For a full discussion of this method in any 
number of dimensions we refer the reader to [H]. For our specific case of one- 
dimensional systems, the phase space density is a function of two phase space 
variables p(x,p) that gives the number of particles dN inside a small phase 
space area da: dp: 

dTV = p{x,p) dx dp . 

In terms of a (large) number of individual particles, this can be though as the 
smoothed-out version of the phase space point density 

N 

p( x ,p) =^2 5 ( x - Xi)5{p-pi) (1) 
i=i 

where (xi,pi) are the phase space positions of the particles. 

We shall consider a collection of non-interacting spinless fermions with a 
single-particle hamiltonian H sp (x,p). For simplicity, we assume that we have 
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identified Darboux coordinates in phase space, with single-particle Poisson brack- 
ets 



Ordinary particle position and momentum are such coordinates, but in principle 
we can use any other canonical pair. 

Semiclassically there is one quantum state per area 2irh in phase space. The 
exclusion principle mandates that we can put at most one particle per state, 
and implies a maximal value for the fermion density p = 1/2ttH, corresponding 
to a dense collection of fermions each occupying one quantum of phase space 
area. 

The ground state of the system would consist of fermions occupying all 
available single-particle states up to a maximal energy Ep. Semiclassically, 
this corresponds to a density p — p Q for all points in phase space with energy 
E = H sp (x,p) < E-p and p = for points with E — H sp (x,p) > Ep. This is the 
'phase space droplet' picture. This state is parametrized by its one-dimensional 
boundary, in this case the set of phase space points at energy Ep: 



We shall assume that this is a single connected curve, although situations where 
the boundary is disconnected, corresponding to topologically nontrivial or dis- 
connected droplets do arise and can be easily treated. The density is given 
as 

p(x,p) = p a d{Ep - H ap (x,p}) , 

where d(y) = i(l + sgn(y)) is the step function. The Fermi level Ep and the 
total number of particles N are related as 



where Ap is the phase space area of the ground state droplet. 

Excitations of the system can be visualized as deformations of the droplet, 
in which its boundary deviates from the Fermi level Ep. The density function 
remains constant and equal to p inside the droplet and zero outside. Such 
states are determined by the (arbitrary) shape of the boundary. A convenient 
way to parametrize the boundary is to give the value of one of the phase space 
coordinates, say p, on the boundary as a function of the other: 



We shall call the function R(x) the boundary field. The density of the droplet 
is given in terms of R(x) as 



The phase space density description of the system has a hamiltonian struc- 
ture, inherited from the single-particle canonical structure. The density function 



{x,p} sp = 1 . 



{(x b ,p b ) : H sp (x bl p b ) = E F } . 




p b = R(x) . 



p{x,p) = p 'd{R{x) -p) . 



(2) 
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satisfies the field theory Poisson brackets 

{p(x,p),p(x',p)} FT = d x p(x,p)S(x - x')6' (p - p') - d p p(x,p)6' (x - x')6(p - p') . 

(3) 

These are the standard Poisson field brackets, implying the canonical structure 
for test functions p[f] = J dx dp f(x,p) p(x,p): 

{p[f], pMh-r = pl{f,g}s P ] 

The hamiltonian is the total energy H of the state, found by summing the 
energies of all particles: 

H = J H sp (x,p) p(x,p)dxdp = p[H sp ] (4) 

The above reproduce the evolution of the field p as the canonical equation of 
motion 

P = {p,H}ft = {p,H sp } sp 

This is the same as the equation implied by the evolution of the particles making 
up the state according to their single-particle equations of motion. 

For a droplet state, the density inside the droplet remains constant, by Li- 
ouville's theorem. Its evolution is determined by the evolution of the boundary, 
that is, the boundary field R(x). The hamiltonian becomes 

H = p J H sp {x,p)d(R(x)-p)dxdp (5) 

The canonical structure for the field R{x) can be determined through a hamil- 
tonian reduction of the original Poisson brackets © to the submanifold of func- 
tions p(x,p) of the droplet form. The result is 

{R(x), R(x')} FT = -—S'(x - x') = -2nh8'(x - x') 

Po 

Conversely, the above Poisson structure for R(x) implies the Poisson brackets 
(j2J for the density function p(x,p) as expressed in j2]). 

The advantage of the above semiclassical construction is that its main fea- 
tures survive quantization. Specifically, the boundary field's Poisson brackets 
are elevated to quantum commutators (from now on we use the convention h = 1 
and denote quantum operators with hatted symbols): 

[R(x), R(x')} = -27ri S'(x - x) (6) 

We recognize the above commutation relations as a chiral abelian current alge- 
bra; the field R(x) is the chiral field that bosonizes the fermion theory around 
the Fermi level Ep. If the fermion system has more than one Fermi points, 
the values of the boundary field at each point become independent commuting 
chiral fields. 
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The above chiral algebra has a unique irreducible unitary representation. 
Its states are into one-to-one correspondence with the excitations of the fermion 
system around the Fermi level. Furthermore, the hamiltonian JSj, diagonalized 
in the states of this representation, reproduces the correct excitation energies 
of the fermion system, at least for the cases of free relativistic or nonrelativistic 
fermions. We recover the results of standard (abelian) bosonization, with the 
boundary field R(x) becoming the chiral part of the many-body particle current 
operator in the bosonic description. For free fermions on a compact space 
there are generically two Fermi points and the corresponding boundary fields 
constitute the chiral components of the derivative of the boson field. 

The mapping of states fails only when the number of excited fermions reaches 
N, in which case we deplete the Fermi sea. This is a large-deformation effect in 
the droplet picture, expressed by a condition of the form R(x) > 0, which would 
introduce nontrivial constraints in the states of the Hilbcrt space of R(x) and 
would complicate the picture. Such effects are purely nonperturbatively in the 
number of fermions N and can be ignored for large N. In the sequel we shall 
talk about the bosonization scheme being 'exact' and will always mean up to 
nonperturbative effects. 

3 The algebra of density operators 

The quantum density operator p{x,p) would naively satisfy the commutation 
relations implied by elevating the Poisson brackets into quantum commuta- 
tors. This, however, is not the case. To see this, start with the quantum version 
of the many-body density expression : 

N 

p( x >P) =^2 S ( X - Xi)8(jp - pi) . 

z=l 

Since Xi and pi do not commute, the above expression is not even hermitian 
and requires the specification of a quantum ordering. 

It is convenient to adopt the Weyl (symmetric) ordering between Xi and pi. 
In terms of the Fourier transform p{l, k) of p(x,p) this is expressed as 

N 

i=i 

Using the Baker-Campbell-Hausdorff formula, it is now a simple matter to 
derive the commutators 

nw — i'k\ 

[p(l, k),p(l', k')} = 2isin I J p(l + l',k + k') . (7) 

This is the famous 'sine' algebra [S] that deforms the classical field theory Pois- 
son brackets Q. It is realized, for instance, by the spatial density of planar 
particles in the lowest Landau level (modulo a factor arising from the different, 
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anti-normal ordering in that case) In a different basis, it is known as the 
Woo algebra JU|- We shall simply call it the (quantum) density algebra. 

The use of the above Weyl-ordered operator p(x,p) is that, if the single- 
particle quantum hamiltonian H sp is expressed in a Weyl-ordered form in terms 
of x and p (which is always possible), then the exact quantum many-body 
hamiltonian can be expressed as 

N . f 
H =}^H ap (xi,pi) = / H sp (x,p)p(x,p)dxdp, 

i—l 

with H sp (x,p) the Weyl 'symbol' of H sp , that is, the Weyl-ordered classical 
expression of H sp in terms of x and p. In effect, the many-body hamiltonian is 
obtained via the (properly ordered) classical expression by the use of p(x,p). All 
other many-body operators can similarly be obtained as appropriate expressions 
involving p{x,p), which becomes a 'universal' many-body operator. 

In the context of the phase space bosonization, what remains to be shown is 
that the droplet expression @ for the quantum density operator p(x,p) in terms 
of R(x) satisfies the proper quantum mechanical commutation relations. If this 
is the case, then the bosonization of any many-body fermion system has been 
achieved. Indeed, as argued before, R(x) reproduces the proper many-body 
fermion Hubert space and p{x,p), acting on this Hubert space and satisfying 
the proper algebra, would be the correct quantum mechanical expression for the 
density operator. In particular, the expression @ for the hamiltonian would be 
exact and would reproduce the correct energy states. 

It is possible to start with the expression (J2J in terms of R(x) and, using the 
operator commutation relations 10 for R(x), show that it formally obeys the 
correct operator algebra |12| . This proof, however, suffers from shortcomings 
related to the presence of singularities. Specifically, the expression for p(x,p) 
is formal and contains infinities. A proper definition requires normal ordering 
and the normal-ordered expression for p(x,p) fails to reproduce the correct 
algebra. 

In the next sections we show that a proper quantum mechanical definition of 
the droplet density p(x,p) requires a nontrivial modification; once this is done, 
the correct operator algebra is recovered. 

4 The semiclassical droplet density operator 

Although the semiclassical droplet density presented above is not the appropri- 
ate expression for the quantum density operator p(x,p), it will nevertheless be 
used as the starting point for an exact bosonization of the system. 

For concreteness, we consider fermions on a periodic space of period L = 2tt, 
which allows for discrete integer Fourier modes for functions of x. We shall use 
Darboux coordinates (x,p) € [— 7r,7r] x [0, oo) where the momentum is positive. 
(This amounts to considering only 'half the Fermi sea states, and thus working 
around one Fermi level with one boundary field operator.) So (x,p) can be 
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thought of as angular and radial variables describing the droplet. The quantum 
boundary field R(x) obeys the commutation relations |JBJ, or, in terms of modes, 

[Rn, Rm] = nSn+m , (8) 

where 

R n = ^- [ R(x) e~ inx dx (n E Z) 
2n J 

denote the Fourier components of the field R(x). 

The above is a set of mutually commuting harmonic oscillators. In standard 
current algebra representation fashion, we consider the Fock space of the oscilla- 
tors consisting on the vacuum state |0), which is annihilated by all the positive 
modes R n (n > 0), and all the excited states generated by the action of creation 
operators R— n . Since the zero mode i?o is a Casimir operator, its action on the 
vacuum -Ro|0) = |0) defines a conserved quantity N, which corresponds to the 
number of fermions. For concreteness, let us assume that N is an odd integer. 
We shall use the notation 

|m, ...,Uj) = R- ni ■ ■ -R- nj \0) 

to represent the action of the negative modes; one should observe that, since 
negative modes commute among themselves, any permutation of (n l7 ... ,rij) 
defines the same vector, so the state \n\, . . . ,nj) is only defined up to permuta- 
tions. 

A first attempt to define the quantum density operator would be to normal 
order the expression for il}(R(x) —p), so as to have an operator with finite matrix 
elements between excited states; that is, 

which is tantamount to moving all the positive (resp. negative) modes to the 
right (resp. left) in any polynomial in the Fourier components of R(x). As the 
positive and negative modes commute among themselves, the above expression 
is meaningful. We can thus introduce the quantum counterpart 

Q(x,p) = ^-:ti(R(x)-p): (9) 

of the semiclassical density. 

It can be seen, however, that this is not the correct quantum expression, 
as it does not satisfy the proper algebra. In fact, it is instructive to study the 
spectrum of the quantum hamiltonian operator 

H e = J H sp (x,p)g(x,p)dxdp-E . (10) 

and compare with the exact many-body fermion spectrum. Here the integral is 
taken over the domain [— 7r, n] x [0, oo) and we have subtracted the inessential 
central term 

E = — / H sp (x,p) i3(R -p) dxdp 
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to shift to zero the ground state energy. Let us henceforth assume that the 
single particle Hamiltonian is 

H sp (x lP )=p s (11) 

for some positive integer s. With this assumption, one can easily find the 
expression 

^ Rni ■ ■ ■ Rn s (12) 
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for the integral 1)10(1 . Since Rq = one has the series expansion 

iRn 

|s(s - 1)(y) (R- n R^ m R n+m + R 

— n—m 



H e R-nRn 



(13) 

Here and in what follows, all the dummy indices n, m range from 1 to oo unless 
otherwise stated, and summation over repeated indices is understood. 

In order to study the low energy modes of H g , it is convenient to define the 
order of a state \nx, . . . , rij) of the Fock space to be the integer n\ + ■ ■ ■ + rij. 
From Eq. (|12f> it is clear that the subspace of states of any fixed order r is 
invariant under H g . Hence both |0) (ground state) and |1) (first excited state) 
must be eigenstates of H e , whose corresponding energies 

E§ = 0, E^sfi)'- 1 (14) 

can be easily obtained from Eq. i|13[l . 

Let us compare energies (|14|l with the actual ones. The exact energies can 
be obtained by quantizing the single-particle momentum, which is restricted to 
take integer values, and then filling the energy levels with N fermions. Hence, 
the ground state energy is obtained by filling all the momenta up to p = 
with fermions. The first excited level corresponds to lifting the top fermion (at 
p = ^y^) to the next available level, i.e., p = The second excitations are 

obtained by promoting either the top fermion to p — or th e second highest 
(p = ^j^) to p = ^y^- Subtracting the ground state energy, this yields the 
values 

£° xact = , (15a) 
E* aet =2-'[(N + l)'-(N-iy], (15b) 

^exact = 2 - S [±{N± 3) S + {NT I)"] ■ (15c) 
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for the lowest eigenenergies of the system. Therefore, the Hamiltonian IjlUfl 
solely provides the leading scmiclassical term of the exact energy. 

We therefore conclude that the normal-ordered semiclassical density d(R(x) — 
p) is, as expected, exact only in the semiclassical limit N —> oc. A further mod- 
ification is needed for the exact operator. 

5 The exact quantum density operator 

The main result of this paper is the construction of a modified density operator 
p(x,p) which satisfies the sine algebra 

ny — i'k\ 

[p(l, k), p(l', k')] = 2isin (-2- 1 ) Kl + l',k + fc') (16) 

and provides an exact bosonization of the system via 

H = J H sp (x,p) p(x,p)dxdp. (17) 

As before, p(l, fc) denotes the Fourier transform of the density operator, in terms 
of which one can recover p(x,p) as 

p{x, p) = ^ E eilx I fc ) e ~ ikp dk ( 18 ) 

/— — oo 

We claim that the modified density operator is given by 

p(/, fc) = / l —r e a{x ' k) e Hx ' k) e ikko e~ ilx dx , (19) 

J 4:TTi sin | 

where 

a(x, k) = — sin — R^ n e~ inx , b(x, fc) = - sin — R n e inx . (20) 
n 2 n 2 

Since the classical limit amounts to fc, I <C 1, and the Fourier transform of the 
step function is -4, it is clear that we recover g(x,p) in the classical limit from 
Eq. (HHJl. 

Let us prove that the modified density operator satisfies the sine algebra . 
We shall multiply each term of the Fourier series in Eq. 1)19(1 by e~2 n to render 
the sum (strongly) convergent, and take the limit e I at the end of the cal- 
culation. This is equivalent to acting with the operators p{l, fc) on the vacuum, 
or states generated from the vacuum through the action of a finite number of 
creation operators R- n , which is the proper physical Hilbert space. On such 
states only a finite number of terms in the expansion of p(l, fc) survive and all 
expressions are finite and convergent. To simplify the notation, from now on we 
shall write a, a' instead of a(x, fc), a(x' ', fc') (and similarly for b). 
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Since the commutator 



■' a 2 Sm g_ „(-ix-f)» (is'-f)ro 



[a, b'} = -4nS nm 2 2_ e ( - ix -5)"e 

n m 

'(1 - e i ( !L ^-+ x - x '+^))(i - e '( 

log ' 



log 1 + 



(1 _ e K h ^T-+x-x'+ie)\(i _ ^(-S-fS.+x-x'+ie)\ 

2 sin £ sin k ' 



2 2 

cos _ cos(x — x' + ie) 

is a central element, we can use the Baker— Campbell— Hausdorff formula to write 

r ~\(k+k r )RQ—\(lx+l'x') 

[p(l,k),p(l',k')} = - / [e"e 6 ,e"e 6 ] — — - r — g dzdx' 
7 167T 2 sm | sm y 

e a+a' e fc+fc' f( x _ x >^ fc ' ; /') e i(fc+fc')Ao-i(*+i> da; (Ja/ ; 

(21) 

where the function / = f(x — x' , k, k' , I') is given by 

Jil'Vx— x') 

f = 1 E p-fe^-e^) 

16tt 2 sin | sin V y 

e il'(x-x'} ( \ 1 



cos cos(a; — x' — ie) cos — cos(x — x' + ie) 

To evaluate the commutator <|21|) , it is convenient to change variables in the 
integral from x and x' to x and y = x — x' , and use periodicity to evaluate the 
y-integral over the boundary of the rectangle 

n = {y e C : -7T < Rcy - ^f- < n, (sgnZ') Jmy > 0} . 

For concreteness, let us suppose that I' is positive. Let us denote the operator- 
valued function e a+a e b+b by F and set K = ^4^-. The function fF has two 
poles in f2, namely y± = ±K + ie, and its residues are 

e il 'y± 

Res(/F;y±) - ± 8^h77?^ 1 -- (22) 

In the limit e J, 0, the poles tend to ±if and one can write 

pi | _ +k _ e f fl_„G- i " x (sm^+sin^o ±i " K ) e f-i? m c imx (sin^+sin^c Tim/f ) 



a(x^~2" ,k+k') b(x^f ,k+k') 



(23) 
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Combining Eqs. I|22[l and 12MI) one can evaluate the integral H21[) by residues as 



[p(l,k),p(l',k')] =27rilim / [Res{fF;y+) + Res{fF;y-)]c i ( k+k ' )A °- i V +l '*> x dx 



This completes the proof of Eq. J7J. 

6 The energy spectrum 

We shall now demonstrate that the spectrum of the Hamiltonian i|17|) expressed 
in terms of the above density operator encodes the full low-energy dynamics of 
the iV-fermion system. 

By an appropriate canonical transformation, almost any single-particle hamil- 
tonian can be brought to the form 



where x is periodic and p > 0. This is true for hamiltonians with compact, 
simply laced equal-energy lines {-ff S p = E}, where x is the angular variable 
around equal-energy lines and p an appropriate conjugate variable. Corre- 
sponding quantum hamiltonians can also be brought to this form upon proper 
reordering. 

It is easy to show that all such bosonized many-body hamiltonians commute 
with each other. Indeed, the expression 



will involve only the Fourier modes /5(0, k). All such modes commute with each 
other, and so will the H h . They constitute a space of commuting 'charges', 
with the hamiltonians H s arising from the single-particle hamiltonian h{p) = p a 
forming a basis for this space. All such hamiltonians have the same set of 
eigenstates, which are, then, independent of h or s. 

We shall, therefore, consider such a hamiltonian H s , and, for simplicity, will 
also subtract the zero-mode contribution to the energy: 




H sp (x,p) = h{p) , 




(24) 




(25) 
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A more manageable formula for H s can be obtained performing the sum over 
I and the p-integral in Eq. I|25|l , yielding 



E 



/— — oo 

oo oo 

E E 

I— — oo i,j=0 

1 d 



167r 3 i sin 



_ e ik(R„-p)-ii(x~x') dx dx i dk dp _ Eq 



pS e ik(%-p)+ily 

87r 2 ii!j! sin | 



a l hP dx dy dk dp - E 



Airi dk s 



fe=0 



e 2 



E 

i,j=0 



1 da; 



(26) 



Despite the presence of a sine in the denominator, it is not difficult to check 
that the integrand is in fact nonsingular at k = 0. 

Now we shall compute the lowest energy levels of this Hamiltonian. From 
Eqs. i|20[l and l|26|l it stems that H preserves the subspace of states of a given 
order because of the a;- integral. (In fact, the order of the state is simply the 
eigenvalue of the Hi operator.) One should also notice that the action of a term 
/ a % W dx over a state of order r vanishes whenever i or j are greater than r. As 
a consequence of these observations, |0) and |1) must be eigenfunctions of H s 
and |0) lies in its kernel. 

The above observation also implies that the sole summand in Eq. I|26|l whose 
action on the order-one state can be nonzero is 



dx ab II) 



-87rsin 2 §|1) 



This easily leads to the formula 



Ei = 2(-iy- 



-l 9 s 
dk s 



k=0 



(sin § e 1 - ) 



^ \1r - 1 



r=l 



N s 



= 2- s [(n + iy -(n- iy] 

for the first excited energy level. This result coincides with the exact en- 
ergy (jl5b|) . 

It can be verified that higher energy levels are also exactly recovered. For 
simplicity we shall concentrate on the second excitations of the system, which 
can be obtained by diagonalizing the action of the Hamiltonian on the linear 
span of |2) and |11). We shall prove that they do coincide with the exact 
energies (|15c|l . 

By the preceding remark, the summands in Eq. I|26|l with i,j > 2 yield 
a vanishing contribution. In fact, a short calculation using the commutation 
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relations © shows that the only terms which contribute to the integral are 

J dx ab 1 2) = -4tt sin 2 k |2) , (27a) 
r a 2h 

/ da;— |2) = -87risin 2 § sinfc |11> , (27b) 

y dxo&|ll) = -167rsin 2 | |11) , (27c) 



dx— — |11) = 167rsin 4 § |11) , (27d) 



4 

ab 2 



dx — |11) = -87risin 2 | sinfc |2) . (27e) 

Using Eqs. I|2t)|l and (|27|l one can write the action of the Hamiltonian on the 

state 1 2) as 



^/isinfc |2) _ 2sin fe Bink \ n) \ 
fc=0 V sinf J 

= 2 ~ S E (o - ?} NS ~ 2r+1 (l + 3 2r_1 ) |2) + 2- s ]T ( * ) J >2 ' 

1 — n V / 1 — n V / 



r=0 v 7 r=0 



|AT s -^(3^ r -l)|ll). 
(28) 

The analogous expression 

L 2" J / \ L 2" J / \ 

H\\\) = 2- s Y^ L _ JiV s - 2r+1 (l+3 2r - 1 ) |ll)+2- s ^ fj7V s - 2r (3 2r -l)|2), 

can be obtained directly from the above equations using the elementary identity 
sin § — sin 3 % — f' n \ . Therefore the action of H is diagonal in the basis 

l 2 4 sm = ° 

^j(|2) ± |11)) and the corresponding eigenvalues are 
£ 2 ,± = 2- s ]T ( r S _\N s - 2r+ Hl + i 2r - 1 )±2- s ^(^)N s - 2r (3 2r -l) 

r=0 ^ ' r=0 \ ' 

= 2- s ^( S )iV s -'-[±(±3rT(Tin 

= 2~ s [ ± (TV ± 3) s T (-/V T ■ (29) 

Thus the exact result (|15c|l is reobtained, as we wanted to prove. 

Notice that the eigenstates corresponding to the eigenvalues ()29fl just ob- 
tained are, up to a normalizing factor, |2) ± |11). The coefficients are not 
incidental: they can, in fact, be read off from the irreducible representations 



13 



of the symmetric group. Higher excited states of order n can similarly be con- 
structed as characters of representations with n boxes in their Young tableaux. 
A complete proof of this statement is possible but will not be given here. 

The simplest examples of physical systems of the form are those with 
linear (s = 1) and quadratic (s — 2) dispersion relations. The s = 1 system 
can be interpreted as a harmonic oscillator in the action-angle basis, or as a 
massless relativistic particle on a circle. The s = 2 system would describe a free 
nonrclativistic particle on a circle. 

Interestingly, Eq. i|14f) shows that the semiclassical approximation is exact 
for these systems. This is related to the fact that the expression for H s obtained 
from 1251) is of the form 

H s = -y / : R(x) s+1 : dx + lower order terms (30) 

(s + 1). J 

where the extra terms include powers : R(x) s ~ 1 : or lower, as well as terms with 
derivatives of R(x). Such terms are all trivial for s = 1, 2 and therefore the 
semiclassical result is essentially exact. 

7 Conclusions 

We have derived the droplet expression for the fermion phase space density 
satisfying the proper operator algebra. We also verified that it reproduces the 
correct many-body fermion energies for non-interacting hamiltonians and that 
it leads to the standard relativistic boson field theory and cubic collective field 
theory for relativistic and nonrelativistic free fermions respectively. 

The main thrust of this work, of course, is the higher dimensional situa- 
tion. For this, the corresponding expression for the higher dimensional phase 
space density should be obtained. Again, a naive (normal ordered) extension of 
the corresponding semiclassical expression is not appropriate and a nontrivial 
quantum modification would have to be performed. 

This problem is presumably related, although not equivalent in any obvious 
way, to the problem of identifying the Fermi operator. In ordinary (abelian) 
one-dimensional bosonization this is given by the exponential of the boundary 
droplet field. This expression is not directly motivated by any droplet consid- 
erations, and does not generalize to the nonabelian case. Similarly, its general- 
ization to higher dimensions is an open issue. This and other questions are the 
subject of further investigation. 
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